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Abstract 



o 

> 

The present article considers time symmetric initial data sets for the vacuum Einstein 
\ ' field equations which in a neighbourhood of infinity have the same massless part as that 

^ . of some static initial data set. It is shown that the solutions to the regular finite initial 

£f>l • value problem at spatial infinity for this class of initial data sets extend smoothly through 

the critical sets where null infinity touches spatial infinity if and only if the initial data sets 
coincide with static data in a neighbourhood of infinity. This result highlights the special 
role played by static data among the class of initial data sets for the Einstein field equations 
whose development gives rise to a spacetime with a smooth conformal compactification at 



5— I , null infinity. 
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q ■ 1 Introduction 

O 

\Q . The analysis of the structure of spatial infinity using the Einstein conformal field equations and 

^sO ' the construction known as the cylinder at spatial infinity — see [8j [TTJ [TH HH [20l [211 [22] — suggests 

that static initial data sets play a special role among the class of initial data sets for the vacuum 
Einstein field equations whose development has a smooth conformal compactification at null 
infinity. This issue is of fundamental importance in the construction of so-called asymptotically 
simple spacetimes. In connection with this expectation, the following theorem has been proved 
in [22]: 



Theorem 1. Consider a time symmetric initial data set for the Einstein vacuum field equations 
which is conformally flat near infinity. The solutions to the regular finite initial value problem 
at spatial infinity is smooth through the critical sets where null infinity touches spatial infinity if 
and only if the data is exactly Schwarzschildean in a neighbourhood of infinity. 

The context of this theorem is better appreciated if one recalls that the Schwarzschild space- 
time is the only static spacetime with conformally fiat slices — see [TT]. The use of conformally 
flat data sets in the analysis of [22] is a convenient technical assumption. Otherwise, the full 
complexity of time symmetric initial data sets would make it impossible to identify useful struc- 
tures. The analysis in [22] builds upon the original analysis in [8] and the computer algebra 
calculations of [19j , to make generic assertions about a certain type of asymptotic expansions for 
the conformal Einstein field equations made possible by the framework of the cylinder at spatial 
infinity. 

The purpose of the present article is to consider an extension of Theorem [T] to the case of 
more general — non-conformally flat — time symmetric initial data sets. Again, explicit computer 
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algebra calculations carried out in [3D] suggest the type of result that one can expect to prove. 
A fundamental difficulty in the way of constructing a generalisation to Theorem [1] is to obtain a 
paramctrisation of time symmetric initial data sets for which it is simple to decide whether the 
data under consideration is static or not. 

The properties of time symmetric initial data sets for the vacuum Einstein field equations, 
(S,hij), to be analysed in the present article are best discussed using a point compactification 
of the initial hypersurface S and an associated conformally rescaled metric hij — the conformal 
metric. The 3-manifold S contains singled out points i 2 ,...} representing the points at 
infinity of the 3-manifold S. Our analysis will be local to one of these points, so it will be 
assumed, without loss of generality, that there is only one of them. This point will be deonted 
by i. In [SJ it has been shown that static initial data sets satisfy a certain regularity condition 
involving the Cotton tensor and its higher order derivatives — see equation ([6]) in the main text. 
This property is, however, not enough to fully assert whether an initial data set is static — initial 
data sets like those of Misner [T3] and Brill-Lindquist [3] satisfy the regularity condition as they 
are conformally flat, but clearly they do not, in general, give rise to static spacetimes. Further 
conditions need to be imposed on the initial data to obtain a static development. This gap 
between initial data sets satisfying the regularity condition and initial data sets which are exactly 
static at spatial infinity is still to be understood^. 

The present article gets around the difficulty exposed in the previous paragraph by considering 
a restricted class of initial data sets for which it is simple to decide whether they are actually 
static or not. This class of initial data sets is constructed by looking at solutions to the equation 
giving rise to the conformal factor $ relating the 3-metric hij and hij, the Yamabe equation — 
see equation (j4]) in the main text. In a suitably small neighbourhood of infinity, the solutions 
to the Yamabe equation can be split into its massless and massive parts. The massless part 
contains the information of the local geometry in a neighbourhood of i, whereas the massive part 
contains information of global nature — in particular the mass. The class of initial data sets to 
be used in the present article takes the solution of the Yamabe equation for static data and adds 
to it a further massive term that does not contribute to the mass so to obtain a new solution 
to the Yamabe equation — this can be done because of the linearity of the setting. This new 
solution to the Yamabe equation implies, in turn, a new solution to the constraint equations in 
a neighbourhood of infinity with the same conformal metric hij as a static initial data set. We 
say that these solutions to the constraint equations have a static massless part. It can be verified 
that this class of initial data sets satisfies the regularity condition of [6] — cfr. also equation (O 
of the main text. This observation is of relevance in the present article as it has been shown in 
[5] that this condition is necessary for solutions to the regular finite initial value problem of the 
conformal Einstein field equations to extend smoothly through the critical sets where null infinity 
touches spatial infinity. 

For the class of time symmetric initial data sets for the Einstein vacuum field equations 
discussed in the previous paragraph one can prove the following generalisation of Theorem [TJ 

Theorem. Given an initial data set with static massless part, the solution to the regular finite 
initial value problem at spatial infinity is smooth through the critical sets if and only if the data 
is exactly static in a neighbourhood of infinity. 

In other words, the smoothness of the development through the critical sets forces the extra 
massive part added to the (background) static initial data to vanish. A more precise version of 
this result will be given in the main text. 

As in the case of the assumption of conformal flatness made in [32], the use of initial data 
sets with a static massless part is a useful technical assumption which allows to identify relevant 
structures in the conformal field equations. It is clear that not all time symmetric initial data 
sets admitting an analytic conformal compactification at infinity have a static massless part. A 
general version of the main theorem of this article can only be obtained once one knows what 

Recently, there has been some progress in this direction — H. Friedrich, parallel session talk in the GR19 
Conference, Mexico. 
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extra conditions have to be imposed on a generic time symmetric initial data set to have a static 
massless part — cfr. similar remarks in the previous paragraphs. This task requires learning 
how to exploit to the maximum extent the conformal gauge freedom implicit in the conformal 
Ansatz. As the conformal metric, h^, encodes all the freely specifiable information of a time 
symmetric initial data set, the extra requirements will have to be in the form of conformally 
invariant conditions on the conformal class. 

The proof of the main theorem builds upon the analysis of the conformally flat case discussed 
in [35]. This analysis relied heavily on the use of computer algebra methods to transform the 
transport equations implied by the conformal field equations at the cylinder at spatial infinity into 
a system of reduced ordinary differential equations for which explicit solutions can be computed 
for any order of the expansion. The approach in the present article consists in conveniently 
grouping the various terms appearing in the transport equations in Schwarzschildean terms and 
deviations-from-Schwarzschild terms. The former are formally identical to terms appearing in 
[2"2"] and thus, assertions about their smoothness can be readily given. It turns out that most 
of the terms that one needs to consider are Schwarzschildean terms. It is only in the last step 
of the argument that deviation terms arise. As it will be seen, their presence indicates that the 
extra massive part that has been added to the static data has to vanish up to a certain order 
— thus, putting into action an inductive argument from which the main theorem is obtained. 
Remarkably, essentially all the computer algebra required for this argument has already been 
performed in [22] , 

Outline of the article 

Section [2J discusses some general properties of time symmetric solutions to the Einstein constraint 
equations in the conformal setting. It also introduces the class of time symmetric initial data 
sets that will be used in our subsequent analysis — initial data sets with a static massless part. 
Section [3] gives a concise summary of the framework of the cylinder at spatial infinity and of the 
so-called F-gauge. Its main purpose is to introduce the notation to be used in the rest of the 
article. It also provides an overview of the key properties of the transport equations implied by 
the conformal Einstein field equations at spatial infinity. Section 0] briefly discusses the key result 
of the construction of the cylinder at spatial infinity for static spacetimes — namely, that the 
structures are as smooth as they can be. Section [3] discusses key properties of initial data sets 
which are static up to a certain order — the results will be used extensively in the sequel. Section 
[5] contains the core of our analysis: a discussion of the properties of solutions to the transport 
equations at the cylinder at spatial infinity for data which is static up to a certain order. The 
results presented in this section take the form of an induction argument which leads, ultimately, 
to our main theorem in Section [7] 

Notation and conventions 

The present article is concerned with the asymptotic properties of spacetimes (M,g^) solving 
the Einstein vacuum field equations 

R»v = 0. (1) 

The metric g^ will be assumed to have signature (+, — , — , — ) and /i, v, . . . are spacetime indices 
taking the values 0, . . . , 3. The spacetime (M. ) g flv ) will be thought of as the development of a 
time symmetric initial data set (S,hij) where S is an asymptotically Euclidean hypersurface. 
The metric hij will be taken to have signature (—,—,—). The indices i, j, . . . will be spatial ones 
taking the values 1, 2, 3. The spinorial conventions of [THUS] will be adopted. The present article 
draws heavily from the analysis in [HI [HJ H2] so we have followed the notation and conventions 
of these references as closely as possible. 
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2 A class of time symmetric data 



For time symmetric initial data sets (<S, hij) the Einstein vacuum field equations ([1} imply the 
constraint equation 

f = 0, on S, (2) 

where r denotes the Ricci scalar of the metric hij . 

Our analysis will be local to a neighbourhood of infinity. Hence, without loss of generality, only 
one asymptotically flat end will be assumed. The asymptotic flatness of the time symmetric initial 
data (S, hij) will be expressed in terms of conditions on a conformally rescaled manifold. For this, 
it will be assumed that there is a 3-dimensional, orientable, smooth compact manifold S, a metric 
h{j, a point i e S, a diffeomorphism $ : S \ {i} —> S and a function G C 2 (S) n C ca (S \ {i}) 
with the properties 

n(i) = 0, D j n(i)=0, DjD k n(i) = -2h jk (i), (3a) 
0>Oon5\{i}, (3b) 
hij=n 2 $*hij, (3c) 

where Dj denotes the Levi-Civita covariant derivative of the 3- metric hij. For the sake of sim- 
plicity the last condition will be written as hij — £l 2 hij so that S \ {i} is identified with S. If 
assumptions (|3a[) - (|3cp are satisfied, the pair (S,hij) will be said to be asymptotically Euclidean 
and regular. Suitable punctured neighbourhoods of the point i are mapped to the asymptotic 
end of S. 



2.1 Asymptotically Euclidean and regular data 

The Hamiltonian constraint, equation ([2]), together with the boundary conditions (|3aJ) - (|3cp imply 
on B a (i) the Yamabe Equation 

h - V) i? = -47r<y(i), d = n~ 2 , (4) 



where S(i) denotes the Dirac delta distribution with support on i while A and r correspond, 
respectively, to the Laplacian and the Ricci scalar of the conformal metric hij . For later use we 
define 

LU = ■ ■ 5 

y/\D k QD k n\ 

It is well known — see e.g. [5J [TT] — that if a is suitably small, then the solutions to the 
Yamabe equation Q on B a (i) are of the form 

d = + W, \X\ = ((X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 ) 1 /^ 

Fl 

for some asymptotically Cartesian coordinates x % . The terms f//|x| and W will be known, re- 
spectively, as the massless and massive parts of The function U, the Green function, satisfies 
the equation 

\ h - \r\ (£) = -47r«(<), 



8 J \\x\. 

and describes the local geometry in B a (i). The function W satisfies the equation 

A h -lr\W = 0, 



annd contains information of global nature. In particular, W{i) = m/2, where m is the ADM 
mass of the initial data set. 
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A rescaling 

h,, ■ > i,',, w '/-,,. n^n' = e 2 n, 

with a smooth positive factor 9 satisfying 8(i) — 1 leaves the physical metric hij = £l~ 2 hij 
unchanged. However, it implies the transitions 

It' I 

<& ^ = O^tf, U h> U' = r-rO^U, W ^ W = B~ l W, 

\x\ 

where \x'\ is given in terms of the /i'-normal coordinates. There are several possibilities to fix this 
conformal gauge freedom. For the purpose of the present analysis it turns out that a good choice 
is that of the so-called conformal normal (cn)-gauge introduced in [8]. In the following definition 
let denote the Schouten tensor of the metric hij . 

Definition 1. The metric hij is said to be in the cn-gauge if given a solution (x(X),b(X)) to the 
3- dimensional conformal geodesic equations 

x u D 1/ x" = -2{b v x v ) x" + (h Xp x x x») W u b Ul 
x»D y b^ = {Kx v ) b p - ±(h x Pb x b p ) h^x v + l Xtl x\ 

x{0) = i, h(x,x)(i) = -l, 6(0) =0, 
(b,x) = 0. 

Remark. Let hij be analytic in a neighbourhood of i. Assuming that a is sufficiently small, there 
exists onS a (i) a unique analytic rescaling h%j i-> — 4 hij for which h\j is in the cn-gauge. The 
metric and connection remain unchanged at i. As it will be discussed in the sequel, the practical 
advantage of the cn-gauge is that it renders simpler multipolar expansions for various quantities 
of interest. 



with initial conditions 
one has that 



2.2 Asymptotically static data 

Let (S, hij) denote a static initial data set given in the cn-gauge, and let 

denote the corresponding solution to the Yamabe equation (U) in a suitably small neighbourhood 
B a (i). The static initial data set can be specified entirely in terms of its multipole moments — see 
e.g. [12]. This important fact will, however, not be exploited here. As a consequence of the 
analysis in [TTJ [T2] one has the following: 

Proposition 1. Let (S,hij) be a static initial data set and let B a (i) be a suitably small neigh- 
bourhood of i. If the conformal metric hij satisfies the cn-gauge B a {i), then there exist normal 
coordinates x % such that hij, U and W are analytic in the neighbourhood. 

One also has that — cfr. [5J [5] — : 

Proposition 2. The Cotton-Bach tensor bij of the conformal metric hij of a static initial data 
set (S,hij) satisfies the regularity condition 

C(Dfc,---Z> fcl 6 w )(i)=0, q = 0, 1,.... (6) 

Remark 1. Given a sequence of multipoles for a static solution, one can always assume without 
loss of generality, that the dipolar terms vanish — this amounts to the gauge choice of working 
in the centre of mass. It can be readily verified that this assumption has the consequence that 
static data in the cn-gauge satisfies 



5 



The latter form of the function W will be assumed in the sequel. 

Remark 2. If U — 1 and W = to/2 in B a (i), then the static initial data set corresponds to 
initial data for the Schwarzschild spacetime. 

The following observation will be crucial in our subsequent analysis. Let W satisfy 

Aft - ~rW = 0, W(i) = 0. 
Clearly, due to linearity one has that 

{) = ^- + W + W (7) 

M 

is also a solution to the Yamabe equation Q with the same boundary conditions in B a {i). Due 
to the analyticity of hij, the function W will also be analytic. Time symmetric initial data sets, 
(iS, hij), constructed out of the conformal metric hij and the conformal factor i9 given by equation 
([7]) will be said to have a static massless part. Moreover, following the ideas in |22j . we make the 
following ad hoc definition: 

Definition 2. A time symmetric initial data set in the cn-gauge will be said to be static up to 
order p m if and only if there exists coordinates x l in a suitable neighbourhood B a {i) of infinity 
such that the solution "d to the Yamabe equation ((4]) is of the form 

U °° 
#=-r- + W + W, W= J2 ■■■x t » ) (8) 

' ' p=p.+i 

where U /\x\ and W are, respectively, the massless and massive parts of a static initial data set. 

Remark. It can be shown that in the cn-gauge the polynomials 

are harmonic with respect to the flat Laplacian. That is, the constant tensors 

w il---i P ,+n w ii---i P ,+2> w ii---i P ,+3 

are trace- free with respect to the fiat metric Sij . 

2.3 Properties of data which is static up to a certain order 

Consistent with Definition [2l we write 

= + Cl, uj = lj + LU, 

where Cl and w are, respectively, the static parts of O and uj (obtained by setting W = 0). Their 
non-static parts Cl and Co are then obtained via 

= — 0, LO = LU — tjj. 

The following result is obtained by direct computation observing expression ([7]): 
Lemma 1. Let (S,hij) be an initial data set which is static up to order p m > 2. Then 

Cl = 0{\x\ p '+ 4 ), lo = 0{\x\P'+ 3 ). 

Let bij denote the Cotton-Bach tensor of the conformal metric h^ of an initial data set which 
is static up to a certain order in the sense of Definition [51 Crucially, one has that: 
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Lemma 2. Let (S,hij) be an initial data set which is static up to order p m > 2. Then 

C(D lq ■ ■ ■ DiJ)ij)(i) = 0, g = 0,... 
where C denotes the operation of taking the symmetric trace-free part. 

Proof. As already discussed, Proposition [5] shows that the property holds for exactly static initial 
data sets. In [8] it has been shown that this property concerns only the massless part of time 
symmetric data sets. The result follows then by observing that an initial data set which is static 
up to a certain order has the same massless part as the one of a static initial data set so that 
bij = bij. □ 

3 The cylinder at spatial infinity 

In [8] a representation of the region of spacetime close to null infinity and spatial infinity has 
been introduced — see also the comprehensive discussion in [TT] . The standard representation of 
this region of spacetime depicts i° as a point. In contrast, the representation introduced in [8 
depicts spatial infinity as a cylinder — the cylinder at spatial infinity. This construction is briefly 
reviewed for the case of time symmetric initial data sets with an analytic conformal metric in a 
neighbourhood B a (i) of infinity. The reader is referred to [3 [11] for a thorough discussion of the 
details — see also [22] . 

3.1 The Manifold C a 

Starting from the initial hypersurface <S, the construction introduced in [H] makes use of a blow- 
up of the point i £ S to the 2-sphere § 2 . This blow-up requires the introduction of a particular 
bundle of spin- frames over B a - Consider the (conformally rescaled) spacetime [M^g^) obtained 
as the development of the time symmetric initial data set (S, hij). Let SL(S) be the set of spin 
dyads 8 = {Sa}a=o,i 011 & which are normalised with respect to the alternating spinor €ab in 
such a way that eoi = 1- Let r = \/2eo, where eo is the future g-unit normal of S and Taa' 
its spinorial counterpart. The spinor taa' enables the introduction of space-spinors — sometimes 
also called SU(2) spinors, see [TJI2JQ2]. ^ defines a sub-bundle SU(S) of SL(S) with structure 
group SU(2,C) and projection w onto S. Given a spinorial dyad S £ SU(S) one can define an 
associated vector frame e a , a = 1,2, 3. We shall restrict our attention to dyads related to frames 
{ej}j = o.... t s on B a such that e$ is tangent to the /i-geodesics starting at i. Let H denote the 
horizontal vector field on SU(S) projecting to the radial vector e^. The fibre 7r _1 (i) C SU(S) 
(the fibre "over" i) can be parametrised by choosing a fixed dyad 8* and then letting the group 
5/7(2, C) act on it. Let (—a, a) 3/jh> S(p 7 t A B ) £ SU(S) be the integral curve to the vector H 
satisfying 5(0, t A B ) = 5(t A B ) £ 7r _1 (z). With this notation one defines the set 

C a = {S( P ,t A B ) £ SU(B a ) | \ P \ < a, t A B e SU(2,C)}, 

which is a smooth submanifold of SU(S) diffeomorphic to (—a, a) x SU(2,C). It follows that the 
projection map 7r of the bundle SU(S) maps C a into B a . The manifold C a inherits a number of 
structures from B a . In particular, the solder and connection forms can be pulled back to smooth 
1-fornis on C a satisfying the structure equations which relate them to the curvature form. The 
explicit form of the structure equations will not be required here. As S is 3-dimensional, the full 
content of its curvature can be expressed in terms of the spinorial counterpart of the trace-free 
Ricci tensor sabcd and the Ricci scalar r. 

In the sequel t A B £ SU(2, C) and p £ R will be used as coordinates on C a . Consequently, one 
has that H = d p . Vector fields X±, X relative to the SU(2, C)-dependent part of the coordinates 
can be introduced by requiring the commutation relations 

[X,X+] = 2X+, [X,X_] = -2X_, [X + ,X_] = -X, 

and by requiring that they commute with H = d p . More importantly, it can be seen that for 
p £ B a \ {i} the projections of the fields d p , X± span the tangent space at p. Given these vector 
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fields, define the frame cab = C {AB) by 

cab = c\ B d p + c\ B X + + c^X-, 

where 

i + _ 1 .+ _ _ 1 

C AB — X AB, Cab ~ ~p ZAB ' C AB' C AB ~ ~p^ AB C AB' 

with constant spinors x ab, Dab and zab given by 

1 i i 1 

XAB = V' Ze (A e B) ' V AB = ?K e A e B > Z AB = —7E e A~ e 



AT 1 _ 1 1 x 

: V 2e (A 6 B) ' ^-B = _ "^7| 6 A e B ) Z AB = ~/=e 

and analytic spinor fields satisfying 

For the frame cab, the connection coefficients, jabcd, are given by 

"lABCD = -1*ABCD + 1ABCD: 1*ABCD = ^ACXBD + ^BDXAc), 
P z 

where 

Iabcd = 0{p), JnCD = 0. 

3.2 An orthonormal basis for functions on SU(2,C) 
Given t A B e Sf7(2,C), define 

/ \ 1/2 / \ 1/2 

T 1 (t A ) - ( m \ ( m \ t (Bl ■■■t Bmh 
^o°o(^s) = 1' 

with j, fc = 0,...,m and m = 1,2,3,.... The subindex expression (A!---A m ) k means that the 
indices are symmetrised and then k of them are set equal to 1, while the remaining ones are 
set to 0. Details about the properties of these functions can be found in [SJ [S] . The functions 
s/m + IT \ form a complete orthonormal set in the Hilbert space L 2 (p, SU(2, C)), where p 
denotes the normalised Haar measure on SU(2,C). The action of the differential operators X± 
on the functions T m • is given by 



X + T m * - ^j{m- ] + l)T m *_ l5 X_T m * = - y/(j + l)(m-j)T m k j+1 . 
In the sequel, we will need to linearise products of the form T m \. To this end, we use the formula: 

A" 

T i 1 3 \ 1 x T i 2 n k 2 = ^2 D ( i i,ji,k 1 ;i 2 ,j 2 ,k 2 ;i 1 + i 2 - 2p, ji + j 2 -p,&i + &2 - p) 

xT J'i+ia-p cq\ 

with /Lt = mm{ii, i 2 , ji + 32, ki + k 2 } and 



D(ix,ji, kr, i 2 , j2, k 2 ;l, m, n) = C(ix,3i\h,fe; I, m)C(i 1 ,ki;i 2 , k 2 ;l, m), 
and C(i,j; k, I; m, n) the Clebsch-Gordan coefficients of SU(2, C). 
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3.3 Normal expansions at X' 



In the sequel, we will consider the lift of analytic fields defined on B a to C a - In particular, the 
lift of \x\ is p. More generally, let £a 1 b 1 ---a 1 b 1 denote a spinorial field on B a . Denote, again, by 
^AiBi— AjBj hs lift to C a - Denote by £j = £,(a 1 b 1 ---a 1 Bi) - ^ < j < I its essential components. The 
function £j has spin weight s — I — j and a unique expansion of the form 

oo p+l 2q 

£j = £j,pP P > G,P = C?' ; p;2g,fe^2o 2q-/+7> 

p=0 g=max{|Z— p} fc=0 

with complex coefficients Cj,p;2g,fc- More generally, we shall consider symmetric spinorial fields 
(,A!—A r onC a with independent components £j = ^(^ 1 ...A 2r ) j , < j < 2r, and spin-weight s = r— j 
which do not descend to analytic spinor fields on B a . In this case one has that 

oo q(p) 2q 

= ^^^3,pP t 0,P = £j,P\2q,kT 2a "a-r+i- 

p=0 q=\r—j\ k=0 

where one has a priori that < \r — j\ < q(p). An expansion of the latter form will be said to be 
of type q(p). 

3.3.1 Particular results concerning expansion types 

We recall the following result in [8]: 

Lemma 3. Assuming the cn-gauge and if for some non-negative integer k it holds 

d (e p f p ■ ■ ■ D El F 1 b A BCD){i) = 0, p = 0, 1,- - • ,k, 
then we have the following expansion types: 

type(r) = type(W) = p, 

type(sABCD) = P, 

type(7 abc d) =p-l, 

tyP e ( a As) = P~ 1 > 
type(C/ - 1) = p - 2. 

In addition, we have that 

U = l + 0(p 4 ), S A BCD = 0(p 2 ), c± B = G(p 3 ), J A BCD = 0(p 3 ). 

Remark. The particular structure of the expansions dictated by this lemma greatly simplifies 
the subsequent analysis. 

Important for our subsequent analysis is the particular form of the lift to C a of the function 
W appearing in the conformal factor d for data which is static up to order p # . Using the methods 
of [5] one finds that 

p.+3 2p oo p 2q 

E E/^ r 2 //+ E EE/^ ? y. (io) 

p=p.+l k=0 p=p.+4q=0k=0 

where the coefficients w p -2 q ,k are related to the constant tensors Wi 1 ...i p via contractions with 
Infcld-van der Waerden symbols. In particular, it can be seen that 

w P :2q,k = 0, q = 0, ...p, k = 0,...2q^^w il ...i p = 0. 

The key observation in equation (JTDJl is that the terms of order 0(p p,+1 ), 0(p p ' +2 ), 0( / o p,+3 ) 
are formally identical to the ones appearing in the corresponding expansions for the functionVT 
of data which is Schwarzschildean up to order p, (W = m/2). In other words, 



Wp.+l;q,k 


= o, 


q 


= o,. 


. .p. - 


i-i, 


k 


= o,. 


..q, 


Wp.+2:q,k 


-o, 


q 


= o,. 


. .p. - 


1-2, 


k 




..q, 


Wp.+3 ]q ,k 


-o, 


q 


= o,. 


. .p. - 


1-3, 


k 


= o,. 


..q. 
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3.4 The spacetime Friedrich gauge 

The formulation of the initial value problem near spatial infinity presented in [Sj employs gauge 
conditions based on timelike conformal geodesies. The conformal geodesies are curves which 
are autoparallel with respect to a Weyl connection — i.e. a torsion-free connection which is not 
necessarily the Levi-Civita connection of a metric. An analysis of Weyl connections in the context 
of the conformal field equations has been given in [7] . In terms of this gauge based on conformal 
geodesies — which shall be called the Friedrich gauge or F-gauge for short — the conformal factor 
of the spacetime can be determined explicitly in terms of the initial data for the Einstein vacuum 
equations. Hence, provided that the congruence of conformal geodesies and the fields describing 
the gravitational field extend in a regular manner to null infinity, one has complete control over 
the location of null infinity. This can be ensured by making B a suitably small. In addition, the 
F-gauge renders a particularly simple representation of the propagation equations. Using this 
framework, the singular initial value problem at spatial infinity can be reformulated into another 
problem where null infinity is represented by an explicitly known hypersurface and where the data 
are regular at spacelike infinity. The construction of the bundle manifold C a and the blowing up 
of the point i e B a to the set 1° C C a , briefly described in section IXT1 are the first steps in the 
construction of this regular setting. The next step is to introduce a rescaling of the frame bundle 
so that fields that are singular at 1° become regular. 

Following the discussion of [S] assume that given the development of data prescribed on B a , 
the timelike spinor t aa introduced in section I5TT1 is tangent to a congruence of timelike conformal 
geodesies which are orthogonal to B a . The canonical conformal factor rendered by this congruence 
of conformal geodesies is given in terms of an affine parameter r of the conformal geodesies by 

e = K- 1 n(i-'?i], (in 



with uj as given by equation ((5]) and where = d" 2 and d solves the Yamabe equation Q — see 
[5J 110) . The function n > expresses the remaining conformal freedom in the construction. It 
will be taken to be of the form k — n' p, with k! analytic, n'(i) = 1. Associated to the conformal 
factor there is a 1-form d^ from which the Weyl connection can be obtained. In spinorial terms, 
one has that for conformally flat data 

1_ ao b j j _ fUx AB -pD AB U-p 2 D AB W 



dAA ' = V5 TAA ' dT ® ~ T A ' dAB ' dAB = 2p v (u + pwy ) • (12) 

The function k in the conformal factor 0, induces a scaling S A ^ ^^ 2 ^>a of the spin frame. 
Accordingly, one considers the bundle manifold C a>K = K}/ 2 C a of scaled spinor frames. Using C a K 
one defines the set 

M a , K = i(T, q )\ q eC a , K ,-^<T<^ q) 



K{q) n{q) 

which, assuming that the congruence of null geodesies and the relevant fields extend adequately, 
can be identified with the development of B a up to null infinity — that is, the region of spacetime 
near null and spatial infinity. In addition, one defines the sets: 

!={{r,q) £M a , K | p(q) = 0, \t\ < l}, 
X ± = {(r,q) S M a , K | p(q) = 0, r = ±l}, 

^ = j(r,<z) € M a , K | p(q) > 0, r = 

which will be referred to as, respectively, the cylinder at spatial infinity, the critical sets and 
future and past null infinity. In order to coordinatise the hypersurfaces of constant parameter r, 
one extends the coordinates (p, t A B ) off C QiK by requiring them to be constant along the conformal 
geodesies — i.e. one has a system of conformal Gaussian coordinates. 
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Remark. For the purpose of the analysis carried out in this article it turns out that the most 
convenient choice of the function k in the conformal factor of equation (llip is 

k = p. 

This leads to considerable simplifications in all the relevant expressions. From this point onwards, 
this choice will always be assumed. 

3.5 The evolution equations 

On the manifold Ai a .n it is possible to introduce a calculus based on the derivatives d T and d p 
and on the operators X+, X^ and X. The operators d p , X + , X- and X originally defined on 
C a can be suitably extended to the rest of the manifold by requiring them to commute with the 
vector field d T . In order to derive the propagation equations, a frame caa' and the associated spin 
connection coefficients Taa'BC of the Weyl connection V will be used. The gravitational field is, 
in addition, described by the spinorial counterparts of the Schouten tensor of the Weyl connection, 
Qaa'bb 1 , and of the rescaled Weyl tensor, 4>abcd — see [7J|H1[TT]. Let fa = fa\ABCD)t- in the 
present gauge, the information of the spacetime spinors c AA >, Qaa'BB' and Taa'bc is encoded, 
respectively, in space spinors c AB , Oabcd and Tabcd — see [TT] for the detailed relation between 
the two sets of spinors. 

In what follows, we will arrange the independent components of the spinorial fields c AB , 
Qabcd and Tabcd in the vector v, and those of 4>abcd in the vector <fi. Suitable field equations 
for the fields contained in v and <p can be obtained from the first and second Cartan structure 
equations, the Bianchi identity of the conformally rescaled spacetime ( M , g pv ) , and the Bianchi 
identity of the physical spacetime (M.,g pv ) — see e.g. [11] for details. A solution to the 
equations thus constructed implies a solution to the vacuum Einstein field equations. The vector 
v has 45 independent complex components, while cf> has 5 independent complex components. 
Using the F-gauge it can be shown that the extended conformal field equations given in [5] imply 
the following evolution equations for the unknowns v 

d T v = K-v + Q{v,v)+L-fa (14) 

where K and Q denote, respectively, linear and quadratic constant matrix- valued functions with 
constant entries, and L is a linear matrix-valued function with coefficients depending on the 
coordinates and such that L| p= o = 0. For the unknowns <p, the Bianchi identity \7 AA fa\BCD = 
implies, respectively, a set of propagation and constraint equations of the form: 

E-9 T + A^(c)^0 = F(r)-0, (15a) 
B"(c) • = G(r), (15b) 

where E denotes the 5x5 unit matrix, A M (c) and B M (c), /i = 0, . . . , 3, are, respectively, 5x5 and 
3x5 constant matrix- valued linear functions depending on the frame field coefficients c AB . Finally, 
F(r) and G(r) denote constant matrix-valued linear functions of the connection coefficients 
Tabcd- 

3.6 Initial data for the evolution equations 

For quick reference, we include here the expressions for the initial data for the curvature compo- 
nents appearing in equations ([T4" )) - (|15bp . These expressions can be deduced from the conformal 
constraint equations — see [8]. On C a , K one has that: 

k 2 1 
Qabcd = --qD( AB D CD }Q. + —n 2 rh A BCD, (16a) 

K 3 

4>abcd = { d (abD C d) + Qsabcd) , (16b) 

where, ft denotes the lift to C a<K of the conformal factor obtained from equation (|7J|. The ex- 
pressions for the frame and connection coefficients do not involve the function J7, and thus their 
detailed expressions will not be required here. 
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3.7 The transport equations 



Consider now the system (fT^|) - (ll5aD with data given on C a . K , and let be given a neighbourhood 
W of C atK in SA a , K on which a unique smooth solution of the Cauchy problem exists. From the 
point of view of the propagation equations, the subset WHl is a regular hypersurface. Introduce 
the notation 

« (0) = v\wm, (O) = 4>\ W rxt- 

Due to the property L| p =o = 0, equations (fH| decouple from equations (|15a[) and can be in- 
tegrated on W 01 using the observation that the restriction of the initial data to 1° coincides 
with Minkowski data. The solutions thus obtained extend analytically to the whole of I and in 
particular to the critical sets I . The set X turns out to be a total characteristic of the system 
(Tl4l) - (|15a[) in the sense that the whole system reduces to an interior system on X. Moreover, the 
constraint equations (|15bl) also reduce to an interior system on X. As mentioned before, this 
feature is a consequence of the fact that the unphysical metric g^ determined by a solution to 
the conformal field equations degenerates as p — > 0. 

A crucial structural property is that 

E + A°(c) = diag(l + t, 1, 1, 1, 1 - r) on X, (17) 

so that the matrix E + A°(c) which is positive definite for |r| < 1 degenerates at 1^. Under- 
standing the effects of this degeneracy is the main motivation behind the analysis in the present 
article and that in [22] . 

In order to analyse the effects of the degeneracy of the matrix E + A°(c) we will consider the 
hierarchy of transport equations which is obtained by repeated application of the operator d p to 
equations (|14 p - (|15b[) and then evaluation on X. By this procedure one obtains interior systems 
for the quantities 

v (p) = dV>v\x, </> (p) = d^4>\i. 
The transport equations take the following form for p > 1: 

d T v (p) = K • v {p) (*) (Q(u ( *\u (p - J ' ) ) +LW) • (j,(P-^ , (18a) 

j=o 

e • d T cf>^ + (A") (0) • = f(t<°>) • </> (rt 

+ J2 ( P ) (F(rW) • <j> (p - ]) - (A")W> ■ d^ p ~ j A , (18b) 

(B")<°> ■ = G(r(°)) • + ( P ) (G(rW') ■ <j> ip - ]) - (B")W • d^ p - j A . (18c) 

j= i V^/ 

Note that the non-homogeneous terms in the equations (|18al) - (|18cl) depend on v^ p \ <j>( p 5 for 
< p' < p. Thus, if their values are known, then equations (I18ap - (|18b[) constitutes an interior 
system of linear equations for and (fi( p \ The principal part of these equations is universal, 
in the sense that it is independent of the value of p. If the initial data on C a . K for the system 
(Tl4l) - (|15a[) is analytic — as it is the case in the present analysis — then suitable initial data for the 
interior system (|18a[) - (|18bl) can be obtained by repeated p-differentiation and evaluation on X°. 

The language of jets is natural in the present context. For p = 0, 1, 2, . . . and any sufficiently 
smooth (possibly vector valued) function / defined on Ai a . K , the sets of functions {f^°\ ■ ■ ■ , 
j(p) j. on X will be denoted by [/] and referred to as the jet order p of f on X — and similarly 
with X replaced by X°. If u = (v, <j>) is a solution to the equations (|18aj) . (|18bl) and (|18cl) . we 
refer to Jjf* [ u ] as to the s-jet of u of order p and to the data J^o [u] as to the d-jet of u of order 

p. An s-jet Jj [u] of order p will be called regular on X = X U X + U I~ if the corresponding 
functions extend smoothly to the critical sets X^ . 

We note the following result of ;22 which will be used in the sequel. 
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Lemma 4. If the s-jets [v] and Jj? [<p] have polynomial dependence in t for some p > 1, 

then jff\v] has also polynomial dependence in r. 

3.8 Decomposition in spherical harmonics 

Our analysis requires decomposing the entries of the vector unknowns and 4>^ in terms of 
the functions T,\. Following the discussion in [22], given a vector — (u^, . . . , u^) with 
entries in 1 x § 2 and non-negative integers q and k = 0, . . . , 2q one defines the sector & q ^[u^] 
as the collection of coefficients 

Ui, P ;2q,k = (2<7 + 1) / uf ] T 2 k dn, 

JSU(2) 

where s is the spin- weight of u\ , d/i is the Haar measure of SU(2) and denotes the complex 
conjugate of uf\ Furthermore, one defines 

2q 

6 ? [uW] = |J 6, )fc [uW]. 

k=0 

With this notation, a sector will be said to vanish if & q [u^}. For convenience in the discussion, 
we introduce the following terminology: 

• A coefficient Ui^q^k will be said static if it depends only on the mass m and on the terms 
in the expansions of U and W . We will write Ui tP -2 q .k — u-i,p\2q,k- 

• The coefficient Mi jP; 2g,fe will be said to be Schwarzschildean if it only depends on m. In our 
gauge, Schwarzschildean terms can only arise in the sectors So- 

• A coefficient Ui jP -2 q ,k will be said to be non-static if it depends on the terms in the expansion 
of W. We will write u^p^fc = u^ p . 2q ,k- 

• A non-static coefficient Ui, p -2 q ,k will be said to be a deviation from Schwarzs child if it only 
depends on m and on the coefficients in the leading term of W — i.e. w p .+i ) 2(p.+i).fe, 
fc = Q,...,2(p. + l). 

The structure of the transport equations for the elements of the various sectors & q ^[v^] and 
& q .k[<t>^] has been discussed in detail elsewhere — see [51[TTJ[22]. In particular, the degeneracy 
at the critical sets X ± implied by expression (|17[) is only relevant for sectors with p > 2. Thus, in 
the sequel it will always be assumed that p > 2. If we denote by ' differentiation with respect to 
r, the relevant transport equations are of the form: 

V 'p;2q,k = K ' v p;2q,k + 2Q(t>0;0 : 0, V p . 2q ,k) + h p;2g ,fc, 
(E + A) • 4>' p -2q,k + Ap;2 9 ,fc ' 4>p;2q,k = fp;2q,fc, 
B ' 4>' p -2q,k + Bp;2 9: fc ' <Pp;2q,k = gp;2q,k, 
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where 
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The terms f p; 2g,fc, g P ;2 g ,fc, h p;2(? ,fc are calculated from the left hand sides of equations (|18a[) - (|18c|) 
using formula ©. Detailed expressions for certain values of the multiindex (p; q, k) will be given 
in the sequel. 



4 The cylinder at spatial infinity for static spacetimes 

As mentioned in the introduction, static initial data sets are expected to play a privileged role 
among the class of time symmetric initial data sets with a development which is asymptotically 
simple. This point of view is supported by the following result proved in |llj showing that 
the construction of the cylinder at spatial infinity described in the previous section is for static 
spacetimes as smooth as it is to be expected. 

Theorem 2. For static vacuum solutions which are asymptotically flat, the construction of the 
cylinder at spatial infinity is analytic in the sense that there exists an a > and a frame for 
which all the fields, including the reseated conformal Weyl tensor extend to analytic fields on 
some neighbourhood J\f of I in M Kt a- This statement depends neither on the particular choice of 
the conformal gauge used to prescribe the (analytic) free datum on S nor on the choice of k. 

Essential for our analysis is the following corollary: 

Corollary 1. The solutions to the transport equations (|18a[) - (|18cp for time symmetric static data 
given in the cn-gauge extend analytically through for all orders p. Moreover, the solutions to 
the transport equations are polynomial in r. 

Proof. Once analyticity has been asserted, the polynomial dependence of the solutions with re- 
spect to t follows from an analysis of the possible solutions to the reduced equations — cfr. the 
discussion in section 6.2 and 6.3 of [22]. □ 

5 Properties of data sets which are static up to a given 
order 

In what follows we will discuss some properties of initial data sets which are static up to a certain 
order — in the sense of definition [2] These properties have mainly to do with the multipolc 
structure of quantities associated to these initial data sets. 

Given a quantity u, we will systematically write it as 

u — u + u, 
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where, in the terminology of the previous section, u and u denote, respectively, its static and 
non-static parts. For quantities on the initial hypersurface C a . K we similarly write uq = uq + uq. 

Let as before, c and T denote, respectively, the parts of the vectorial unknown v containing the 
independent components of the frame and connection coefficients. As mentioned in section T3.61 
the expressions for the value of the frame and connection coefficients on the initial hypersurface 
C a ,n are independent of f2. Thus, one has that: 

Lemma 5. For data which is static up to order p, > one has that 

J.p) — !.{p) t(p) _ -hiv) 

L — L i A — A i 

for all < p < p, . 

Crucial for our analysis is the behaviour of the data for the curvature spinors Qabcd and 
4>abcd- It is observed that 

Sl = p 2 -mp a + ^m 2 + <D(p 5 ), (20a) 

D (AB D CD) n = -6me 2 ABCD + 12m 2 e 2 ABCD p 2 + 0{p 3 ), (20b) 

from where it follows that the leading terms in these expressions are Schwarzschildean. The first 
non-Schwarzschildean contributions arise at orders 0(p 5 ) and 0(p 3 ), respectively — see e.g. |13) . 
Using the expansions (I20a|) - (|20b|) together with formulae (|16al) - (|16b[) one obtains after a long but 
straightforward computation the following: 

Lemma 6. For data which is static up to order p, > 2 one has that the curvature fields on C a ^ K 
satisfy 

Qabcd = Qabcd - Qabcd = 0(pP' +2 ), 

<t>ABCD = tpABCD ~ 4>ABCD = 0{p P ' + 1 ). 

Moreover, the terms 

&(P'+ 1 ) o(P.+2) o(P.+3) o(P.+4) 

^ ABCD ' ^ABCDi ^ABCD' ^ ABCDt 

Ifa.+l) I(P«+ 2 ) ^(P'+ 3 ) 

ABCD ' ^ABCD ' ^ABCD ' 

onC a ,K, are of the form of deviations from Schwarzschild data. On the other hand, ^abcd contains 
more general types of deviations. 

An inspection of the explicit computations in (T3] one obtains the following result. 
Lemma 7. For a data which is static up to order p % > 3 one has that 

=*«>), v^=V {2 \ 
L (D =L (D j L (2) = L(2) ; L (3) = f J (3) j 

/iawe orii?/ contributions to the sector ©o- On the other hand, the coefficients of 

t ,(3) = ^0) j 0(2)=0(2) 5 

contain contributions to the sectors ©o a«rf ©2- 

Remark. In the previous lemma L denotes the linear matrix-valued function appearing in 
equation ([T4|. It depends on the conformal factor 9 as given by ([TT]) and on the 1-form dAB as 
given by (TT2"j) . 

A lenghty but straightforward computation using formulae (1111) and (|12[) renders the following 
lemma. 

Lemma 8. for rfaia which is static up to order p m > 3 one /ias i/iai 

Lb.+i) =0 , L(p-+ 2 )=0. 
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6 Solutions to the transport equations for data which is 
static up to a certain order 

In this section we consider a time symmetric initial data set which is static up to order p = p, + 1. 
The data can be parametrised in the form 

^0=^0+ Vq, 4> = 4>0+ 00 

where Vq, 4>o corresponds to exactly static initial data, while Vq — 0, 0o = if W — to/2. A 
similar splitting can be applied to the quantities v^ p \ <// p '. More precisely, 

v (p) =£b) 0(f) = p>0. 

One has the following result. 
Lemma 9. For initial data which is static up to order p = p,, one has that 

v {p) =0, < p < p. + 1, 

and 

0»=O, 0<p<p.. 
In other words, the first deviations from static evolution appear in <fi( p,+1 \ 



Proof. One first notes that 
and 



4 P) = 0. P < P. + 1 



0o P) =O, p<p.. 



One then argues by induction. It is clear that 

&<® = 0, #<°> = 0. 

Now, given p such that 1 < p < p», assume that i)( p ) = and (f)^ p ■* = for < p' < p. 
Substitution of the Ansatze v^' — + 4>^ — (jyj) -\- j = 0, ...,p into the v- 

transport equation (|18ap renders the equation 

5 T u (p) + <9 r £ (p) = K ■ v {p) + K • v {p) 

+ Y1 ( P -) (Q(^ {j) ,^ (p ~ j) )+^ (j) ■ 4> {P ~ J) 

3=0 \- J ' 
3=0 V/ 
3=0 \- J ' 
3=0 

Using the induction hypothesis and 

d T v {p) = K • «W + M (Q(v (j \v {p - j) ) + L U) • (p ^ } ) , 

3=0 ^ 

one finds that 
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As this last equation is homogeneous in the unknown the initial condition Vq = implies 
-jj(p) — o. A similar argument leads to the following equation for 4>^: 

Again, the homogeneity of this last equation, together with the initial condition fyf 1 = renders 
0» = 0. Thus, one has that = 0, <f>^ = 0, < p < p. . The argument can be repeated for 
«(P' +1 > as 4 P ' +1) =0- □ 

The ideas of the proof of Lemma [9] will be used to study with some detail the solutions to the 
sector 6p.+i transport equations for the orders p. + 1, p. + 2, p. + 3 and p. + 4. 

6.1 The transport equations at order p = p, + 1 

From Lemma |9] one knows that = t>(P« +1 ). Thus, at this order one only needs to study 

the solutions to the Bianchi transport equations. 

Substitution of the Ansatze 

V U) ={fU) +ij<J)^ 0(3) = 0« 0<j<p. + 1 

into equations (|18b|) - (|18c[) with p = p. + 1 and using LemmalHlone obtains the following equations 
for 0(P' +1 ): 

y/2E • «9 r (p) + (A^)(°) • d^ p) = F<°) • 
(B")W -<9 M (p) = gK°).0». 

We now focus on the sector © p . + i of these equations. Recalling that (A^)(°), (B^)(°), F^ ', 
only contain the sector ©o, one obtains the matricial equations 

A ' 0p.+X;2(p.+X) + A p. + l;2(p. + l) ' 0p.+X;2(p. + X) = 0, (21a) 
B ' 0p.+X;2(p.+X) + B p.+l:2(p. + l) ' 0p. + l;2(p.+l) = 0, (21b) 

where A, A p . + x;2(p.+x) > B and B p#+1;2 (p.+i) are the matrices given by (|19a[) - ([19bl) . For the sake 
of the simplicity of the presentation, the subindex k has been omitted from these and most of the 
subsequent equations. A lengthy but straightforward computation shows that the initial data for 
these equations is given by 

0p.+i;2(p.+i),fc(O) = (<2o,fe> 2x,fcj ®2,k, ffli.ki ao,k), (22) 
with k = 0, . . . , 2(p. + 1) and 

ao,fe = -y/p,ip. + l)(p. + 2) (p. + 3)ii)p.+i;2(p.+i),fc 1 
ox.fc = -4(p. + 3) a/ (P. + 1)(P. + 2)wp. + i; 2 (p.+i),fc, 

02, fc = -6(p. + 2) (p. + 3)li)p. +1 ;2(p.+X),fc- 

The previous equations lead to the following crucial observation: 

Observation 1. Equations (|21a[) - (|21bp and £/ieir corresponding initial data are formally identical 
to the sector © p .+i [0( p * +1 - ) ] transport equations for data which is Schwarzschildean up to order 
p, + 1; the solutions are, therefore, also formally identical to those obtained in [22]. 

As a consequence of the latter observation one obtains from the analysis in 22 that: 

0p.+X;2(p.+X) = ¥W(T)(1 ~ tY>-\\ + tY'~\ 

with ^p.+i(t) having entries which are polynomials of degree 4 in r. In particular, 

£ P .+X(0) = 0p.+X;2(p. + l)(O), 

as given by (|2"2"j) . 

Combining this analysis with Corollary [T] one finds the following result: 
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Proposition 3. The solution to the transport equations at spatial infinity at order p. + 1 for 
data which is static up to order p, are polynomial in r. Hence, they extend analytically through 

T = ±1. 

6.2 The transport equations at order p = p, + 2 

Using Lemma @] one finds that the solutions i>(p»+ 2 ) to the order p, + 2 ^-transport equations are 
polynomial in r given that and <f>^ for < p < p. + 1 are polynomial in r. However, we 
require more precise information. 

Again, we consider the transport equations (|18al) for the order p. + 2. The substitution of the 
Ansatze 

uW = + v {p) , (p) - (p) + (p) , < p < p. + 2, 

and considerations similar to the ones used for the order p. + 1 lead to the following equations 
f or £(p.+2) and 0(p.+2). 

9 T t;(p-+ 2 ) = K • ^-+ 2 ) + Q(£( \£( p -+ 2 >) + Q(^-+ 2 \£(°>) + L« ■ 4> ( 'P' +1 \ 
V2E ■ d T 4> {p ' +2) + (A")<°> • a Al 0( p '+ 2 ) = F(°) ■ p . +2 + (p. + 2)F« . 0(p.+D 

+ f^-+ 2 ) . - (p. + 2)(A^)( 1 ) • 

(B")(°)a^^+ 2 ) = G(°) • P . + 2 + (p. + 2)G« • <^' +1 ) 

+ G(p.+2) . 0(0) _ ^ + 2 )(B^)W • dp4> {p ' +1 'K 

Again, as a consequence of Lemma [7] one has that 

LW 

contain only contributions to the sector @o so that one obtains directly the following equations 
for the components of & p .+i[v^ p ' +2 \ cf>( p,+2 )]: 

Up.+2;2(p. + l) = K ^p.+2;2(p. + l) + Q (^0;0 , Up. +2;2(p. + 1) ) 

+ Q( w p.+2:2(p. + l), U 0; o) + Li ;0 • 0p.+l;2(p. + l), (23a) 
(E + A) ■ p . +2; 2(p. + l) + A p.+2;2(p. + l) ' 0p.+2;2(p.+l) = F 0;0 ' </>p.+2;2(p. + l) 

+ (p. + 2)Fi ;0 • </>p.+l;2(p. + l) + Fp.+2;2(p.+l) ' 4>0;0 - (p. + 2)Ai ;0 ■ 0p. + l;2(p.+l) , (23b) 
B ' ^p.+2;2(p.+l) + B p.+2;2(p. + l) ' 4>p. +2;2(p. +1) = G 0;0 ' 4>p.+2 + (p. + 2)Gi ;0 ■ 0p,+l;2(p.+l) 

+G p . +2 ;2(p.+i) • <Po-,o - (p. + 2)B 1;0 ■ 0p.+i ; 2(p.+i)- (23c) 

Using Lemma [9j one sees that the components in 

uo ; o, Li ;0 , F 0; o, Go ; o, 

are exactly Minkowskian while the components in 

4>o-,o, Ai ; o, Bi ; o, Fi ; o, Gi ; o, 

are exactly Schwarzschildean — that is, they depend only on the mass m. On the other hand, as 
already seen, the components in 0p.+i ; 2(p.+i) depend only on u>p.+i ; 2(p.+i),fc- This leads to the 
following crucial observation at this order: 

Observation 2. The equations (|23a[) - (123cp are formally identical to the order p, + 2 transport 
equations for initial data sets which are Schwarzschildean up to order p.; similarly, due to Lemma 
[5| the initial data set is formally also of the form of a perturbation of Schwarzschild. 

As a consequence of the previous discussion one can directly use the analysis and results of [22 
to directly conclude that the solutions # p#+2; 2(p.+i), and 0p.+ ; 2(p.+i)> to equations (|23aj) - (j23cj) 
are polynomial in r. Combining this observation with Corollary [T] one has that: 

Proposition 4. The solutions to the order p, + 2 transport equations (I18a[) - (|18cl) are polynomial 
in t and, thus, extend analytically through r = ±1. 
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6.3 The transport equations at order p = p, + 3 

We now adapt the procedure discussed in the previous section to the analysis of the order p # + 3 
transport equations. As in the previous order the polynomial dependence in r follows directly 
from Lemma [31 once one knows that the entries in 0(p«+ 2 ) are polynomial. Further detailed 
information will follow from the analysis of the sector 6 p ,+i- 

The analysis for this order is similar to that for orders p. + 1 and p, + 2. Substitution of the 
Ansatzc 

v (p) =v(p)+ <j(p) , 0<» = 0» + 0» , < p < p. + 3 
into the transport equations (I18ap - (|18cp leads to the following equations for i>(p»+ 3 ) and 0( p * +3 ); 

9r £( P .+3) = K . -(P.+3) + Q(^(0)^b.+3)) + Q(i j(p. +3)^(0)) 

+(p. +3)Q(w (1) ,-o (p - +2) ) + (p. + 3)Q(-o (p - +2) ,- i ) (1) ) + (p. +3)L (1 > • 0( p '+ 2 ) 
+§(p. + 3)(p. + 2)L< 2 > • 0^- +1 ) + I>+ 3 ) • 

V^E0 r #"' +3 > + (A*»)W • /1 #*+ 3 ) = F<°> • p . + 3 + (p. + 3)F( X ) • 0^+ 2 ) 

+i(p. + 3 )(p. + 2 )F< 2 ) • 0^ +1 > + (p. + 3)F^+ 2 > ■ 0« + F^'+ 3 ) • 0^ 

-(p. + 3)(A^)( 1 ) • 3 M #»*+ 2 ) - ~(p. + 3)(p. + 2)(A" ) (2) • d^'+V 
-(p. + 3)(A^)( p -+ 2 ) -d M (1) - (A^)( p - +3) •a Al (o) , 

(B^)(o)^0(p.+3) = 6(0) . p . +3 + (p. + 3 )G« ■ 0^+ 2 ) + i(p. + 3)(p. + 2)G< 2 ) • 0^- +1 ) 
+G (p.+3) . 0(°) + (p. + 3)G<f+ 2 ) • 0« 

-(p. + 3)(B")W • d M 0^+ 2 ) - l(p. + 3)(p. + 2)(B")( 2 ) • d^'+V 
-(p. + 3)(B^)(f-+ 2 ) • - (B^)(f-+ 3 ) • d^°\ 

Using Lemma[7]one obtains the following equations for the components of & p ,+i[v( p ' +3 ' , 0( p,+3 )]: 

^p.+3;2(p.+l) = K ' ^p.+3;2(p. + l) + Q (^0;0 , +3;2(p. + l) ) + Q(^p.+3;2(p. + l) , ^0 ; o) 

+ +3)Q(t>l ; o,Vp. +2 ;2(p. + l)) + (P. +3)Q(t>p. +2 ;2(p.+l),'Ul;0) 

+(p. + 3)Li ;0 • p .+2;2(p.+i) + 5 (p. + 3)(p. + 2)L 2;0 • 0p.+i :2 (p.+i) 

+Lp.+3;2(p. + l) ' 0O;O, (24a) 
(E + A) • 0p. +3;2 (p. + l) + Ap >+3 . 2 (p. + i) • </>p.+3;2(p.+l) = F ;0 • 0p.+3;2(p. + l) + (P. + 3)Fi ;0 • </>p.+2;2(p.+l) 

+ |(p. + 3) (p. + 2)F 2;0 • 0p.+i ;2 (p. + i) + (p. + 3)Fp >+2;2( p >+1) • 0i ;O 

+ Fp. +3;2 (p. + l) • 00:0 

-(p. + 3)Ai ;0 • 0p. +2;2 (p. + i) - (p. + 3)Aj • 0p. +2;2( p. + i) - (p. + 3)A^ • 0p. +2;2 (p.+i) 

-§(p. + 3)(p. + 2)A 2;0 • 0p. + i ;2(p . + i) - | (p. + 3)(p. + 2)A+ • 0p. + i :2 ( p . + i) 

-|(p. + 3)(p. + 2)A2; • 0p.+i ;2( p. + i) - (p. + 3)Ap. +2;2( p. +1) • 0i ;O! (24b) 

^ ' ^p.+3;2(p.+l) + Bp,+3;2(p.+l) 1 </>p.+3;2(p.+l) = G 0; ' <£p. +3;2(p. +1) + (p. + 3)Gi ;0 • 0p.+2;2(p. + l) 

+ \{j>» + 3)(p. + 2)G 2; o0p.+i;2(p.+i) + G p>+3;2 (p. +1 ) • 00:0 + (p. + 3)Gp. +2;2( p. +1 ) • 0i :O 
-(p. + 3)Bi ; o • 0p. +2;2 (p.+i) - (p. + 3)B+ • 0p. +2;2 (p.+i) - (p. + 3)Bj: ■ 0p. +2;2 (p.+i) 
-\{p. + 3)(p. + 2)B 2:0 • 0p. + i ;2 (p.+i) - | (p. + 3)(p. + 2)B^ • 0p. + i; 2 (p.+i) 
-f (P. + 3)(p. + 2)B2: • 0p. + i ;2 (p. + i) - (p. + 3)B p#+2;2( p. +1) • 0' 1;O . (24c) 
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As in the analysis of lower order transport equations, we begin by noticing that the entries of 

v 0;0, Ll;0) F 0; o, G 0: o 
are Minkowskian — i.e. independent of m and W. On the other hand, the entries of 



«1;0, 


^2-Q, 


</>0;0, 


</>l;0, 






Fl;0, 




Al ;0 , 


A ± 


A 2 ;0, 


A 2;0' 






Bl;0, 


B *0' 




R± 

■ D 2;0' 



are Schwarzschildcan in our gauge — i.e. they depend only on m. 

As a consequence of the analysis of the transport equations for the orders p m + 1 and p m + 2 
one has that the entries in 

Up. + l;2(p.+l), 0p. + l;2(p. + l), 4>p.+2;2(p.+l), 

depend only on m and the coefficients $ p .+i ; 2(p.+i);fe- As in the analysis of the orders p, + 1 and 
p. + 2 one has the crucial observation: 

Observation 3. The equations (|24al) - (|24cp and their initial conditions are formally identical to 
the order p, + 3 transport equations for initial data sets which are Schwarzschildean up to order 
P.- 

As in the analysis of the orders p, + 1, p, + 2 and p, + 3, one can use directly the analysis in 
[2"2"] to conclude that the solutions # p .+3;2(p.+i) and p .+3 ; 2( P .+i) to equations (|24a[) - (|24c[) are 
polynomial in r. More generally, one has that: 

Proposition 5. The solutions to the order p, + 3 transport equations (I18ap - (|18cl) are polynomial 
in t and, thus, extend analytically through t — ±1. 

6.4 The transport equations at order p = p, + 4 

Finally, using similar methods, we discuss the order p, +4 transport equations. As in the previous 
order the polynomial dependence in r follows directly from Lemma @] once one knows that the 
entries in 0(p«+ 3 ) are polynomial. The rest of the analysis is more subtle than at lower orders. 
The reason for this is twofold: 

(i) As a consequence of Lemma [5] one finds that $abct) nas a multipolar structure which is 
more complicated than that of data which just deviates from Schwarzschild data. 

(ii) The discussion of the transport equations for the sectors in S p .+i involves non-trivial 
multiplications of the functions T. 3 ,. 

Evidence from explicit calculations suggests that the solutions to the Bianchi transport equations 
at this order are, generically, not smooth at the critical sets. 

Proceeding as in the case of lower order one finds that the substitution of the Ansatze 

v {p) = £(p) + #(?), 0(p) = 0(p) + 0<p<p.+4: 

into the transport equations (|18al) - (|18cl) and taking into account Lemma [8] leads to the following 
equations for i)(P«+ 4 ) and </>( p,+4 ): 

9r £(p.+4) = K . £(p.+4) + Q(£(0) ) £b.+4)' ) + Q(-(p. +4)^(0)) 

+ (p. +4) (Q(i>( 1 ),tJ('" +3 )) + Q(i(' , «+ 3 ) ) u( 1 ))] 

+ I(p. +4)(p.+3) (Q(^ 2 \^- +2 )) + Q(£( p - +2 \£^)) 

+ L(P.+4) . 0(0) + ^ + 4 )Lb.+3) . 0(D 

+(p. + 4)L« • 0(f'+ 3 ) + |(p. + 4) (p. + 3)f/ 2 > • <^-+ 2 > 
+ i(p.+4)(p. + 3)(p. + 2)L( 3 ).^- +1 ), 
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v / 2Ed r (p ' +4) + (A^)(°>d M <^-+ 4) = F(°) • 0p. +4 + (p. + 4)F« • 0< p -+ 3 ) 

+f(p. + 4) (p. + 4)FW • 0(f+ 2 ) + ±(p. + 4)(p. + 3) (p. + 2)F< 3 > • 
+±(p. + 4) (p. + 3)F(f+ 2 ) • 0< 2 ) + (p. + 4)F(f+ 3 ) • 0« + F(p-+ 4 ' • 
-(p. + 4)(A^)( 1 )a^(P'+ 3 ) - i(p. + 4)(p. + 3)(A" )^d^ {l " +2) 
-i(p. + 4)(p. + 3)(p. + 2)(A^)( 3 )a^(P- +1 ) 

-i(p. + 4)(p. + 3)(A^P'+^d^ - (p. + 4)(A")(p-+ 3 >c> a1 0« - (A^)(P-+ 4 )a M 0(°\ 



P.+3) 



+i(p. + 4)(p. + 3)G( 2 ) • 0>-+ 2 ) + |(p. + 4)(p. + 3)(p. + 2)G< 3 > • 

+i(p. + 4)(p. + 3)6(p-+ 2 ) • 0( 2 ) + (p. + 4)6(f+ 3 ) • 0« + 6(f+ 4 ) • 0(°) 

-(p. + 4)(B^)( 1 )9 m 0(p-+ 3 ) - i(p. + 4) (p. + 3)(B^)( 2 )^0(p«+ 2 ) 
-|(p. + 4)(p. + 3)(p. + 2)(B^)( 3 )a M ^- +1 ) 

-|(p. + 4) (p. + 3)(B")(f +2 )d„0( 2 ) - (p. + 4)(B^)(f-+ 3 )a A1 0( 1 ) - (B^P'+^d^°\ 

As in the previous orders, we extract equations for the components of & Pm j r \[v^ Pm+i \ 0( p,+4 )]. 
In this case the analysis is more involved as there are terms involving non-trivial products of 
spherical harmonics. An inspection shows that the terms containing this type of non-trivial 
products are 

F<3) . 0(f+!), f("-+ 2 ) • 0( 2 ), (A^ 3 ) • a M (p ' +1) , (A^)(p-+ 2 ) • 9 M 0( 2 ), 

go) . 0(p-+!), 6(f-+2) • 0(2), (b^)( 3 ) • 9 M (P - +1) , (B")(f+ 2 ) • a M 0(2). 

All these terms contain products of the form T± 31 ki x ^(p.+i)* 2 ^' wmcn usm g formula © can 
be linearised to render 



J 4 fcj x J 2(p. + 1) j 2 



— „ 7^ J1TJ2 , Tl 

L 2(p.+3) ,ii ,i 2 MM ± 2(p.+S) fei+fc 2 L 2(p. + l),ji j 2 ,fci,fc2 J -2(p. + l) fci+fca-1 

, J1+J2-2 
" rc *2(p.-l),ji,i2,fci,fe2- 1 2(p,-l) fci + fe 2 -2' 

with c 2 ( p . + 3), J - 1); 7 2 ,fc 1 ,fe 2) c 2(p.+i)j lj2 ,fc 1 ,fc 2 , Q2(p.-i),ji,j a ,fci,fe some numerical coefficients. Their 
explicit form will not be essential for the subsequent analysis. Thus, the sector S p .+i of, say, 
p(3) . 0(P.+l) i s of the form 

F 3;0 ' 0p. + l;2(p. + l) + C 4,2(|).+l);p.+l ' F 3; 4 ■ 0p. + l ; 2(p. + l) > 

where C^fo.+iv-p.+i denotes a matrix with numerical entries. The term F 3 ; o • 0p.+i ; 2(p.+i) is 
formally identical to the one one would obtain from considering asymptotically Schwarzschildean 
data. A similar analysis can be carried out with the other terms containing non-trivial products. 

The equations for the components of & p ^+i[v^ p ' +4 \ 0( p,+4 )] are: 

^p.+4;2(p. + l) = K ' ^p.+4;2(p. + l) + Q(#0;0, #p.+4;2(p.+l)) + Q(#p.+4;2(p. + l) , ^0 ; o) 
+ (p. +4) (Q( / U 1;0 ,'U p .+3;2(p.+l)) +Q(Up.+3;2(j,.+l),«l;0)) 
+ |(P. +4)(p. + 3) (Q(l>2;0,Up. +2;2( p. + l)) + Q(#j,.+2;20.+l),U2;o)) 

+(p. + 4)Li ;0 • p .+3;2(p.+i) + |(P. + 4)(p. + 3)L 2;0 • 0p. +2;2 (p.+i) 
+|(p. + 4)(p. + 3)(p. + 2)L 3;0 ■ p .+i;2(p.+i) 

+L p . +4;2 (p.+i) • 0o ; o + (p. + 4)Lp. +3; 2(p.+i) • 0i ; o, (25a) 
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(E + A) ■ 0p <+4; 2(p. + l) + A p.+4:2(p. + l) ' 0p.+4;2(p.+l) — F 0;0 • 0p.+4;2(p. + l) + (p. + 4)F 1; Q ■ 0p. +3;2(p. +1) 

+§(p. + 4) (p. + 3)F 2;0 • ^p.+2;2(p.+i) + |(p. + 4)(p. + 3)(p. + 2)F 3;0 • 4.+i;2(p.+l) 
+F p .+4;2(p.+i) • 00:0 + (P. + 4)Fp #+3;2 (p. +1 ) • 0i ;O + I (p. + 4) (p. + 3)Fp. +2;2 (p. + i) • 2;O 

-(P. + 4) ^Al ;0 • 0p. +3;2 (p. + l) + Aj • 0p.+3;2(p. + l) + A[I • 0p.+3;2(p. + l)) 

-f 0. + 4)(P. + 3) (A 2;0 ■ 0p. +2; 2(p.+l) + A Jo • 0p.+2;2(p.+l) + A^.g • 0p. +2;2 (p. + l)) 

-\{V» + 4)(p. + 3)(p. + 2) (a 3;0 • 0p. + i ;2( p. + i) + A+.q ■ 0p. + i ;2( p. + i) + A^Q ■ 0p. + i ;2( p.+i) 

-(p. + 4)A p>+3;2(p . +1) • 0i ;O - |(p. + 4)(p. + 3)A p#+2;2 (p. + l) • 2 ;O 



R 



2(p.+i);p. + 1' 



(25b) 



B ' 0p.+4;2(p.+l) + B p.+4;2(p. + l) ' 0p. +4;2(p. +1) = G 0;0 ' 0p. +4;2(p. +1) + (P. + 4)Gj. ; ' 0p.+3:2(p. + l) 

+ 4)(p. + 3)G 2;0 • 0p. +2;2 (p.+i) + \{p. + 4)(p. + 3)(p. + 2)G 3;0 • p .+i;2( P .+i) 
+Gp. +4;2 (p.+i) • 0o ; o + (p. + 4)Gp. +3 . 2( p. +1 ) • 1;O + |(p. + 4)(p. + 3)Gp. +2;2( p. +1) • </> 2;0 
-(p. + 4) (Bi ;0 • 0p. +3;2 (p.+i) + Bj ■ 0p. +3;2 (p.+i) + B^ ■ 0p. +3 ;2(p.+i)) 
-|(p. + 4)(p. + 3) (b 2;0 • 0p. +2;2( p. + i) + B^o • 0p. +2;2 (p.+i) + B^g • 0p. +2;2 (p.+i)) 

+ 4)(p. + 3)(p. + 2) (B 3;0 • 0p. +1;2( p. +1) + B+ Q • 0p. + i ;2 (p.+i) + B^ • 0p. + l;2(p.+l) 
-(p. + 4)Bp. +3;2( p. +1) • 0i ;O - |(p. + 4)(p. + 3)Bp. +2;2( p. +1) • cf>' 2 . 

+S 2 (p.+i) ; p.+i (25c) 



where 



R-2(p. + l);p. + l = (P. + 4)(p. + 3)C 4l2 (p. + l);p.+l 



g (p. + 2)F 3;4 • 0p.+i ;2 (p.+i) + 2-Fp >+2 . 2( p. +1 ) • 2;4 

|(p. + 2) (A 3; 4 • 0p. + i ;2 (p. + l) + Aat 4 • 0p.+l;2(p. + l) + • 0p.+l;2(p. + l) 

i(p. + 2) (Ap. + 2;2(j,.+l) • 2;4 + A+ +2;2 (p.+l) • 02:4 + A p, +2 ;2(p.+l) ' 02;4 



(26a) 



S 2 (p.+l);p.+l - (P. + 4)(p. + 3)C 4:2 (p >+1 ). p> + 1 • g(p. + 2)G 3;4 • $p.+l;2(p.+l) + ^Gp. +2;2 (p > + 1 ) • 2;4 



-\{V» + 2) (B 3;4 • 0p. + l;2(p.+l) + • 0p. + l;2(p.+l) + Bg. 4 • 0p.+l ;2 (p.- 

+ 2 ) (Bp. +2;2 (p. + l) • 02;4 



p.+2;2(p. + l) 



■ 02:4 + B 



'p.+2;2(p. + l) 



02;. 



(26b) 



The initial data for the transport equations (125b|) - (|25c|) can be written as 

0p.+4,2(p. + l)(O) = 0p.+4,2(p.+l)(O) + 0p.+4,2(p.+l)(O), (27) 

where 0p.+4,2(p.+i) depends solely on m and Wp.+i ; 2(p.+i),fc (deviation from Schwarzschild) while 
0p.+4.2(p.+i) contains contributions from the multipolar structure of the reference static data. 
In order to discuss equations (|25a[) - (l25c[) we introduce a further Ansatz. We write 

^p.+4,2(p. + l) = 6p.+4,2(p. + l) + £p.+4,2(p. + l): (28a) 
0p.+4,2(p. + l) = 0p.+4,2(p.+l) + 0p.+4,2(p.+l), (28b) 

where £> p .+4.2(p.+i) and t>p,+4,2(p.+i) are the solutions to equations (|25ap ~ (125cp with 

R-2(p. + l);p. + l = 0, S 2 (p >+1 ).p >+1 = 0, 0p.+4;2(p. + l)(O) = 0. (29) 

The crucial observation is the following: 
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Observation 4. Equations (|25a[) - (|25cj) and the data ([2T)) satisfying ([2^)) are formally identical 
to the order p, + 4 transport equations analysed in J2Sfj for data which are Schwarzschildean up 
to order p m . These equations have no solution which is smooth at the critical sets. 

From the analysis in [35] it follows that: 

d 1 T 1 $ P . + 4;2( P .+i) = m 3 < +1;2(p . +1) (l - rf'- 15 (l + r)"- 18 

x U p .+4(t) + <p+ +4 (r) ln(l + r) + £~ +4 (r) ln(l - T )j , (30) 

where the entries of tp, ip^ are polynomials of degree 33. Substitution of the Ansatz (|28a[) - (|28b[) 
into equations (]25bp - (l25c[) renders: 

(E + A) • 0p. +4 ;2(p.+l) + Ap, +4; 2( P .+1) ' 0p.+4;2(p. + l) = R-2(p. + l);p. + l, (31a) 
B ' ^p.+4;2(p. + l) + B p.+4;2(p. + l) ' </>p.+4;2(p. + l) = &2(p. + 1) ;p. + 1 , (31b) 

with R2(p.+i) : p.+i and S2(p.+i) ;p .+i as given in (|26a[) - (|26bp . To conclude our argument, one 
needs to analyse the non-polynomial solutions to these equations. The detailed structure of the 
entries in the non-homogeneous terms R.2(p.+i);p.+i and S2( p .+i) :p .+i can be obtained from the 
explicit calculations of [13 . In particular, one has that these terms are polynomials of order 
2p. + 10 in r. More importantly, it can be explicitly verified that 

d 8 T R2(p. + l);p. + l = R 2( p. + 1); p. + 1 (r)(l - tY '- 10 (1 + rf ~ 10 , 
dr S 2 ( P . + i );P . + i = S 2(p . +1);p . +1 (r)(l - r)f- 10 (l + rf -~ w , 

where R- 2 (p.+i); P .+i( r ) and ^>2(p.+i);p.+i(t) have entries which are polynomials of degree 18 in 



r. Thus, the structure of the solutions to equations (I31a[) - (|31b>l) is best analysed if one takes 8 
r-derivatives of the equations. From here arguments similar to those in [22 allow to show that 

C^p.+4;2(p.+l) = ^2;4 Wp. + l;2(p.+l)(l " ^'-"(l + rf^ 11 

x Up.+4(r) + <ft +4 (r) ln(l + t) + 0- +4 (r) ln(l - T )j , (32) 

with <^ p>+ 4(r), <^^ +4 (r) having entries which are polynomials of degree 24 in r. 

Integrating the expressions (l30l) and (|32l) one hnds that the polynomials multiplying the 
ln(l ± r) in p .+4 :2 (p.+i) are of degree 2p. + 14 and have an overall factor of w 3 w p>+1:2 ( )9>+1 )., 
whereas those in <A p .+4 ; 2(p.+i) are of degree 2p, + 10 and have an overall factor W2;4W Pm+1 . 2 ( p ,+i)- 
It follows that the polyhomogeneous terms in p .+4 ; 2(p.+i) and p .+4 ; 2(p.+i) cannot cancel each 
other to produce a </>p.+4 ; 2(p.+i) which is entirely polynomial. Thus, one has the following: 

Proposition 6. The solutions to the order p, +4 transport equations (I25b[) - (l25c[) with data given 
by equation (l2"Tl) have no polynomial solutions unless u>p.+i ; 2(p.+i),fc = 0. If this condition is not 
satisfied the solutions develop logarithmic singularities at r = ±1 and the solutions are of class 

C""(-i,i)nO +3 [-i,iL 



7 The main result 

The discussion in the previous sections is summarised in the following result. 

Proposition 7. Given a time symmetric initial data set which in a neighbourhood B a of infinity 
is static up to order p = p, (in the sense of Definition^), the solutions to the transport equations 
for the orders p = p. + 1, p. + 2, p. + 3 are polynomial in r , and hence, extend smoothly through 
the critical sets . On the other hand, the solutions at order p — p. + 4 contain logarithmic 
singularities which can be avoided if an only if the initial data is, in fact, static up to order 
P = P. + 1 ■ 
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As in the case of the analysis given in |22j . one can use the previous result to implement an 
induction argument in which the explicit computer algebra calculations of [19) play the role of 
the base step to obtain the main result on this article. 

Theorem 3. The solution to the regular finite initial value problem at spatial infinity for time 
symmetric initial data of the form given by Definition [H is smooth through if and only if the 
restriction of the data to TP coincides with the restriction of static data at every order — that is, 
if their jets of order p coincide for all p. Furthermore, the analyticity of the setting implies that 
the data is exactly static in a neighbourhood of infinity. 

As already mentioned in the introduction, in order to complete the analysis of the solutions 
of the transport equations at the critical sets for general analytic time symmetric initial data 
sets, one needs (conformally invariant) conditions which reduce the initial data set to data with a 
static massless part in the sense discussed in Section [3J This problem will be analysed elsewhere. 
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